FEASIBILITY OF USING SATELLITES IN 
NON-EQUATORIAL, ELLIPTIC, 24-HOUR ORBITS 
FOR UNINTERRUPTED COMMUNICATIONS 


by 

VINOD KUMAR JOSHI 



DEPARTMENT OF AERONAUTICAL ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY KANPUR 

JANUARY 1979 



FEASIBILITY OF USING SATELLITES IN 
NON-EQUATORIAL, ELLIPTIC, 24-HOUR ORBITS 
FOR UNINTERRUPTED COMMUNICATIONS 


A Thesis Submitted 

In Partial Fulfilment of the Requirements 
for the Degree of 

MASTER OF TECHNOLOGY 


by 

VINOD KUMAR JOSHI 


to the 

DEPARTMENT OF AERONAUTICAL ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY KANPUR 

JANUARY 1979 



fit ~ 




Ctlf* ; , 


Nu, A 




/> 

0 


U APE 1379 



CERTIFICATES 





* * ± 


O' 

l 



This is to certify that the work 
FEASIBILITY OF USING SATELLITES IN NUN -EQUATORIAL, ELLIPTIC, 
24-HOUR ORBITS FOR UNINTERRUPTED COlDtUNICATIONS has Been 
carried out under my supervision and has not "been submitted 
elsewhere for a degree. 


Jan. 9, 1979* 


\ S' . h 
jy** > ^ 




\ dh e\j/ 


( KRISHNA KUMAR ) 

Assistant Professor 
Department of Aeronautical Engineering 
I.I.T. Kanpur 


d ' U RADU ATK OFFICE 


j 

i 


V, *, ;hesis has berm approved I 
.j :‘v award of I he I..H gree of | 
Oaoiv of Technolo^ v (*vl Tech.) | 
a aac.ould.nce with \ I to I 


r* * u l a t ions of the Indian 
7 ,i u itute of Technology taflpur 


Slued 


■2«> M? 





ACKN OWIE DGEMBNT 


The author wishes to record his gratitude to thesis 
advisor, Dr. Krishna Kumar for suggesting this interesting problem 
and for his constant encouragement and guidance throughout the 
course of this work. 

The author also wishes to convey his thanks to his 
friends for their cooperative help during this work. He also 
wishes to thank Mr. S.S. Pethkar for his excellent typing. Thanks 
are also due to Mr. Warrayam Singh of Computer Centre. 


V. K. JOSHI 



ABSTRACT 


It is proposed to explore the feasibility of using 
24 -hour, elliptic inclined orbits for satellite communications. 

The study has been initiated with an analysis of the effect of 
orbital geometry on the various important aspects related to commu- 
nications satellites, e.g., boosting energy requirements, visibility 
and changes in peak power and total energy required for communica- 
tions. The analysis shows that the use of the proposed non-stationery 
orbit can lead to significant fuel savings. This may however result 
in some increase in communication power requirements. 

The major challenge in the use of the proposed concept 
appears to be in designing a suitable passive attitude control scheme 
which can provide for the changing satellite attitude so as to 
achieve its fixed orientation with respect to the moving line of 
sight. Fortunately, it has been possible to show that the solar 
radiation pressure can be effectively employed for achieving this 
objective to a high degree of pointing accuracy. 
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1 . INTRODUCTION 


'What we are building now is the nervous system of mankind, 
which will link together the whole human race for better or worse, in 
a unity which no earlier age could have imagined... " said Arthur C. 
Clarke who first conceived of communication satellites more than 30 
years ago. 

1 

What is a communication satellite ? The communication satellite 
can be looked upon as a space relay moving in orbit, receiving infor- 
mation from a ground station and transmitting it to another ground 
station. In several respects, space communication is simpler than 
conventional communication which is severely handicapped due to the 
problems like tropospheric scatter, ionospheric radio or under water 
acoustic communications. As a matter of fact, today - 30 years after 
the prophecy by Clarke - the satellites appear to have already become 
an important link in global communication network. 

The various satellite systems which have been proposed for 
communications are i 

(i) Passive reflector satellites, 

(ii) Active satellites in low and medium altitude orbits, and 

(iii) Active satellites in a 24 -hour synchronous orbits. 

The merits and demerits of the various satellite communication 

2-4 

system have been examined by several investigators . Of the Several 
concepts proposed for space communication, the one most commonly used 
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and now almost universally accepted requires placing the satellites 
into perfect geostationary synchronous orbits. The single over-riding 
consideration that seems to have weighed heavily in its favour is 
that satellites that appear absolutely stationary in the sky from the 
earth allow the use of simpler, less expensive ground stations in 
communications network. However, there are some other important 
considerations which appear to have been ignored. This may perhaps 
be attributed to the undue emphasis being laid on the aspect of 
operational simplicity. The technological problems associated with 
the non-st at ionary satellite alternatives and the likely difficulties 
of excessive demand on altitude control in ecliptic orbits may have 
been othor possible factors. 

The synchronous geostationary satellites although ideally 
suited to uninterrupted communications are difficult to achieve and 
maintain particularly over the extended periods of operations. This 
difficulty arises due to several reasons : 

(i) Large orbital deviations are caused due to the various 

perturbing forces. The frequent orbital corrections are 
therefore needed leading to rather stringent station-keeping 
requirements if the full potentialities of the geostationary 
satellites are to be realized. 

(ii) To inject the satellites into the geostationary orbit, 

large multistage rockets with high boosting capability are 
needed. The use of a non-equatorial launching site 
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( unavoidable if the launching has to he done from within India) leads 
to further increase in expenditure of energy, the additional energy 
being required for affecting the necessary final transfer maneuver for 
orienting the satellite orbit to the equatorial plane. This further 
affects the overall economics of the geostationary satellite adversely. 

5 

It may be pointed out that a recent study by Dial and Cooley 
has established that the net on-station weight capability of satellite can 
be increased considerably by using a 24-hour, inclined, elliptic orbit. 
This analysis is of considerable significance as it suggests the 
possibility of affecting major fuel saving through the use of inclined 
elliptic but 24-hour orbit for satellite communication. 

Here, it is proposed to explore the alternative possibility of 
using elliptic, non-equatorial, 24-hour satellites for uninterrupted 
communications. It is felt that this may in general lead to conside- 
rable saving in terms of booster power requirements without compromi- 
sing to any appreciable extent the communication performance chara- 
cteristics. 

It is felt that the apparently simpler and economical proposition 

if technically feasible would open up a multitude of possibilities 

which are particularly significant in the light of growing problem 

6 

related to the overcrowding of the limited equatorial space with 

7 

satellites and of course the associated political/legal implications. 

Hence, it is proposed to study some important aspects related to the 
use of 24-hour inclined, elliptic orbits for satellite communications. 
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Chapter 2 examines effect of orbital geometry on the requirements of 
■booster capability, communication power and energy, and visibility of 
the satellites from ground station. 

It may, however, be pointed out that one of the most challeng- 
ing problems in the use of the proposed orbits for communications is 
likely to be the control of attitude motion due to excitation induced 
by eccentricity. This problem is further complicated by continual 
periodic drift of the satellites towards or away from local vertical 
through the ground station. This brings us to the important problem 
of control of satellite orientation along the line of sight i.e. the 
line joining the satellite and the ground station executing periodic 
oscillations about the local vertical. This phase of investigation 
aims at exploring the possibility of using solar radiation pressure 
for achieving this variable attitude control. The problem has been 
studied in order of increasing complexity. 

In Chapter an attempt has been made to synthesize a suitable 
but simple control law which can enable the proposed solar controller 
to orient the satellite along the local vertical with a high degree 
of pointing accuracy. 

Chapter 4 deals with the more complicated problem of developing 
a suitable, simple solar controller model for achieving the variable 
attitude control of satellite along the moving line of sight. 



2. EFFECT OF ORBITAL GEOMETRY OH BKERG-Y AND 
COIMUITICATIOHS HEQ.UIHEMH1MTR . 


The various aspects of satellite motion likely to affect the 
acceptability of its orbit for communication are 

(i) boost- energy requirement s, 

(ii) visibility of satellite from ground station, and 

(iii) communication power and energy requirements. 

The study has therefore been initiated with an analysis of the 
effect of orbital geometry on these factors. 

2.1. Boost-Energy Requirements ; 

The analysis compares the effect of orbital eccentricity and 
inclination on boost-energy requirements for the various 24-hour 
orbits. Since injection of satellite into the final desired orbit 
is in general preceded by establishing the circular parking orbit 
first, attention has been focused on computing the total velocity 
increment required in transfer of satellite from the parking orbit 
to final desired orbit. Since the objective of analysis is limited 
to only a preliminary estimate of impulse requirement and their 
comparison, the 2 -impulse, minimum energy i.e. Hohmann transfer 
trajectory has been assumed for transfer to the perfect geosynchronous 
orbit. Furthermore, the final impulse at the apogee of the transfer 
orbit is modified to provide for the change in inclination of the 



6 

orbital plane. In the alternative case of transfer to the general 
24-hour elliptic orbit, a single impulse is assumed to act at the 
perigee (Fig. 2.l) of the desired elliptic orbit. 

The dimensionless expressions for the total velocity- 

increment required in transfer from the parking orbit to the final 

circular or elliptic orbits with common 24-hour period are found 
* 

to be j 


A7 C - 2n 2 “ 2n 2 /( 1+ n 2 ) “ ^2 



Using these expressions, it is possible to compute the percentage 
saving involved in the use of the proposed elliptic orbits instead of 
the perfect circular orbit. These results are shown in Table 1. 

It is thus apparent from these results that the use of elliptic 
inclined orbits can enable a significant saving in booster power 
requirements. The saving increases with increase in eccentricity as 
well as inclination. 


* 


See the Appendix 



PERCENTAGE SAVING IN PO\/ER 


c 

■H 


4 ' 


U J 


C 


cr 
« — . 




rc 

4*' 


KM 

t>- 

4“ 

4" 

CD 



CM 

LTV 

O 

vo 

KM 

Cv 

0 

* 

• 

• 

« 

♦ 

f 

KM 

L T\ 

O 

vo 

T- 




ON 

CTn 

00 

00 


>v 














r •* 







4m 







*H 

4 .. 


CM 

4" 

tr- 

KM 

KM 

4 * 


KM 

CC 

UM 

LTV 

t— 


UM 

• 

• 

• 

m 

O 

*4 

CM 

4- 

00 

km 

00 

KM 



om 

CO 

00 

r>- 

c — 

is 







0 







*H 









\>- 

t — 

CM 

KM 

vo 

*% 


OM 

KM 

O 

T— 

vo 

SS 

O 

» 

• 


* 

* 

O 

CM 

CM 

VO 

O 

4“ 

CD 

Hi 


CM 

00 

00 

c— 

VO 

&H 







<4 







iS 







M 














O 

1 

t-— 

LTV 

vo 

4“ 

cc 

is 


0 - 

CM 

4~ 

■'vl- 


H 

UM 

• 

* 

• 

• 

♦ 


T~ 

0 

CM 

4 - 

r— 

T~ 



OM 

CO 

tM- 

VO 

vo 



T— 

LT\ 

00 

LT\ 

vo 


O 

UM 

t— 

vo 

CM 

KM 


T*” 

* 

* 

* 

• 

* 



VO 

4" 

4" 

VO 

Cr\ 



(O' 

c— 

vo 

l C\ 

4- 



00 

CM 

ON 

4" 

KM 



G\ 

VO 

c— 

CM 

00 


LT\ 

• 

• 

4r 

* 

* 



4~ 

MO 

km 

UM 

OM 



1 

UM 

4" 

KM 

CM 



1 

UM 

CM 

T“ 


4" 


O 

4~ 

00 

T— 

KM 

4~ 



* 

* 

• 

» 

• 



CM 


!>- 

om 

T“ 










l n 

0 

LT\ 

0 

um 

CD 


O 


*r~ 

CM 

CM 



• 

• 

* 

• 

• 



O 

0 

O 

0 

O 

• 







O 

h*>» 










CM 

KM 

-4“ 

UM 

to 








CO 

4" 

t— 

0 

x — 


CM 

KM 

UM 

0 

vo 


« 

« 

• 

• 

• 


KM 

CTn 

LfM 

CM 

00 


tM- 

VO 

vo 

VO 

UM 

1 


CTn 

CM 

4" 

UM 

1 

t-— 


N 

CTn 

OM 

CM 

00 


• 

• 

• 

m 

• 


OV 

•'vi- 

O 

r— 

KM 1 


VO 

vo 

VO 

UM 

UM | 


CM 

4- 

T — 

UM 

j 

KM 


VO 

O 

OM 

CM 

O 


« 

• 

• 

» 

« 

O 

km 

CTn 

4* 

T— 

00 

O 

VO 

um 

LTN 

UM 

4- 

T — 

CTn 

vo 

4* 

v- 

T - 

, t> Q 

vo 

ov 

OM 

vo 

OM 

<} 

m 

• 

♦ 

* 

• 


UM 

O 

VO 

KM 

O 

CD 

UM 

UM 

4 - 

4“ 

4- 

*> 

<3 

1 

T — 

CD 

KM 

vo 

00 

4- 


CO' 

vo 



vo 

11 

• 

• 

• 

♦ 

m 

km 

OM 

VO 

4" 

CM 

O 

4“ 

KM 

KM 

KM 

KM 

l}» 

<3 

KM 

KM 

UM 

4* 

CD 


vo 

OM 

CM 

CM 

VO 


* 

* 

♦ 

• 

• 

** 

vo 

4“ 

4“ 

4“ 

4- 

U 

CM 

CM 

CM 

CM 

CM 

O 

g 

4 

0 

OM 

O 

4“ 

Ov 

* 1 — i 

LTN 


-4- 

CM 

OM 


* 

• 

m 

• 

• 

bn 

rO\ 

UM 

tM~ 

OM 

O 

sT 


n — 



CM 

•H 

> 

Cj 

CQ 

0 

LTN 

O 

UM 

O 

O 

KM 

KM 

4“ 

* 

4- 

• 

UM 

% 

O 

O 

O 

O 

O 

-P 

d 

CD 

O 

U 

0) 

Ph 

VO 


00 

CM 

O 

T— 

1 * 



8 


2.2. Visibility of Satellites from Ground Station . 

An important condition for use of satellites as communication 
relay without time delay ia its visibility from the ground stations 
being used for the purpose. This imposes constraints on the maximum 
permissible orbital eccentricities that can bo used with satellites 
in 24-hour orbits. An attempt has, therefore, been made to obtain a 
theoretical estimate of this upper limit constraint on eccentricity. 
By analyzing the geometry of motion (Fig. 2.2) it is easy to see 
that the visibility of satellite from the ground station is assured 
provided the following inequality constraint is satisfied 

cos 0 > r e /r ... (2.3) 

where, cos 6 = r . r 0 /( |r | |r e | ) ... (2.4) 

On evaluating the expressions for cos 0 in terms of the 

8 

orbital elements, this inequality takes the form 

r R cos u cos(ft-(5 ) - sin u cos i sin ( fi- <$)} cos cj> 

+ ( sin u sin i) sin ^ H > 1 ...( 2.5) 

where r n = r/r Q 

XX = (0 4 . 0 

6 = 5 0 + (0 e t 
«o “ 

♦ 


longitude denoting the initial position of the ground 
station with respect to vernal equinox 
latitude of the ground station 



“e* “C- e sin e A-e 2 / (l+e cose) + 2 tan"V(l-e)/(l+e) tane /2} J 

(2*4) 

This constraint can also he written as 

r n Cc° S (u *-“ e t +Q— S Q ) cos <j> + { (l- cos i) sin ( SJ — 6 Q - w e t)sin u} 
cos$+ sin u sin i sin 4Q > 1 ... (2.7) 

It is apparent that the visibility of the satellite ■would, he 
influenced by it.s initial longitudinal position with respect to that 

of the ground station. Normally, the relative position of satellite 

injection into the final orbit would be so chosen as to ensure the best 

visibility conditions, leading to a considerable simplification ill the 

above relation 

r n [p os ( 0 -« e t) c°s $ - (l- cos i) cos <J> - | sin i sin <f> | J> 1 

- * * • (2 *q) 1 

Using the expression for a> e t obtained earlier, it can he shown that 


cos 


( 0 - w G t) Zfc *1 - 4e 2 sin 2 8 > 4e^ 


... (2 .9) 


Substituting (2.9) into (2.8), the inequality takes the 


form 


/ 


1 - 4e > 1 - cos i + | sin i tan <j> | + sec<J> / (a n (l~e)} 


. . .( 2 . 1 #) 


where 


a n “ V 3 


« 6 .628 
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The result is of considerable general significance being 
independent of the several orbital elements, e.g. 0 , oi > ft t <5 0 * 
Furthermore, it provides useful design information giving the maximum 
permissible eccentricities, consistent with visibility, as a function 
of orbital inclination with equatorial plane as well as the latitude 
of the ground station. Results are presented in Table 2. 

From those results, it is apparent that the use of large 
eccentricities is possible for communications even for relatively 
large inclinations and for large ground station latitudes. For 
example, it is possible to ensure visibility for eccentricities upto 
O.35 even in cases where i and <j> become as large as 30°. 


2.3. Communication Power Requirements : 

Recognizing that the communication power required for trans- 
mission from or to satellite is proportional to the square of the 
distance between tho satellite and the ground station, it can be 
shown that 


P 


max 





... ( 2 . 11 ) 


where 


P = Communication Power required for geosynchronous orbit, 
ref 

P = Maxi mum Communication Power required in case of alternative 
max 

elliptic, 24-hour orbit. 

This expression can be rewritten in a more convenient form as s 



MAXIMUM FERI.HSSIBLE ECCENTRIC HIES FOR UNINTERRUPTEI) VISIBILITY 
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W S »f “D+»/(l-l/^)3 2 ...(2.12) 

or 

W P rof - G + 1 -1T6 e I* ... (2.12a) 

The results obtained using this relation are shown in 
Table 3. It may be noted that the communication power increases 
with increase in eccentricity. However, this increase in peak 
power requirement may not be able to offset the gains of fuel 
saving as it would also be accompanied by the corresponding 
decrease in the power required at the points of closest approach. 
A more realistic comparison would, therefore, be provided by 
estimate of the total communications energy required* 

2 .4 ■ Communication Energy Requirements s 

A comparison of energy requirements for communications for 
the various eccentricities can be made using the relations 

2 it 

B “ / P dfi ... (2.13) 

0 

where the power P can be represented by the relation 

P a (r - r e )^ (2*14) 

Substituting (2.1 4) into (2 .13) and carrying out the necessary 
manipulation, we get 
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Ea J (r 2 - 2r r e + r|)de (2.15) 

0 

On evaluating the integrals appearing in this expression, it is easy 
to show that for elliptic orbits 

E e o 2n [a 2 - 2a r e A- e 2 + r Q 2 U (2.16) 

For circular orbits* i -e * when e = 0, this expression reduces to, 

E c a 2* [a 2 - 2a r c + r e 2 2 (2.17) 

On dividing (2.16) "by (2 .17) we get 

E c /E c = + ( 1 -Vi- e 2 )/(a n -l) 2 (2.18) 

or 

E c /E c - V-e 2 + 0.032(1- M-e 2 ) 

Table 4 shows the effect of orbital eccentricity on increase 

\ 

in cominunication energy in elliptic orbits. It may be noted that 
the total energy requirement shows a slight decrease with increase 
in eccentricity. The gains offered by this marginal decrease in 
communication energy are however not very significant. The requirement 
of increased peak power for communication does appear to an additional 
problem. However if the uninterrupted communications is not essential 
there is a possibility of combining the satellite communication 
during the periods when the satellite is relatively closer to the earth 
with other mission objectives. In such a situation, the main 
advantages of elliptic orbits namely the higher peak power require- 
ments for communication can be substantially overcome. 



attitude control op satellites in elliptic orbit s using 

SOLAR RADIATION" PRESSURE 


5.1 Introduction 

For many space applications such, as communications, weather- 
forecasting, survey of earth-resources, it is necessary to maintain 
the satellites at a fixed orientation with respect to the earth. 

The use of active attitude control methods although desirable for 
achieving high degree of pointing accuracy leads to penalty in terms 
of increased weight, reduced space- availability on board and limited 
satellite life-time . That perhaps explains why there have been 
numerous attempts at developing attitude control methods utilizing 
environmental forces . 

The possibility of using Solar Radiation Pressure (SRP) for 

9-14 

attitude control has been explored by several investigators . Even 
though effectiveness of this concept has been clearly established, 
particularly for satellites in high altitude orbits, it appears to 

have found only limited acceptance by the spacecraft designers. This 
lack of enthusiasm for the solar attitude controllers may be partly 
attributed to the availability of the well proven conventional passive 
and semi-passive alternatives using gravity gradient, satellite spin, 
magnetic torquing etc. The complexity of the proposed controller models 
end their operation has perhaps been a major detrimental factor as 
well. It seems that from the very beginning, there has been undue 
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emphasis on using the combination of proportional and derivative 
feed- hack control policies for achieving high degree *f pointing 
accuracy. In the process, certain obvious possibilities emphasizing 
the simplicity of control and operation appear to have been overlooked. 

This chapter investigates the feasibility of developing a 
simple solar attitude controller for satellites operating in elliptic 
orbits. Attention is focused on evolving a simple control strategy 
which can counter the adverse effect of eccentricity normally 
responsible for worsening of the attitude characteristics of the 
conventional methods. 

3.2 Formulation of the Problem : 

The formulation begins with the analysis of motion of 
a cylindrical satellite in elliptic orbit, in the ecliptic plane 
about the oarth-conter 0 (Fig. 3 . 1 ). The satellite model considered 
here carries two sets of plates made of light, rigid but highly 
reflective material (o.g. aluminized mylar membrane). The two sets 
are assumed to be mounted with their surfaces normal to the orbital 
plane. Denoted by x,y,z are the principal coordinate axes of the 
satellite with origin at the mass-centre s. The angle made by the 
z-axis with the local vertical defines librational angle, ip . Using 

Lagrangian formulation, the governing equations of the planar libra- 

J _. . ... 10,11 

tions can bo written as 

(1 + e cos 6) ip" - 2 e(l + ^')sin e+ 3 % sin $ cos 4> =(n 5 /vl)Q 

... ( 3 . 1 ) 
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where represents the generalized force due to the SEP. 

For specular reflection assumed and using the principle of virtual 

work, the expression obtained for the generalized force takes the 
form 

S " ( s /c' ) B 1 +p 1 ) il 1 A- 1 sin( e + Tp _ <j>)| sin (6 +ip _<f> ) | 

- 0 + V 2) ^2^2 COs ( 0 + ^ ~ $) I cos (0 +1|) _ 4 > ) I 3 

... ( 3 . 2 ) 

substituting the above expression for in (3.1) transforms the 
equation of motion to 

(l + e cos 0 ) 4 ’" - 2e (l + t|>') sin0+ 3 ]L sin^cos’J' = 

■z * ' 

Ui + o)7(i + e cos 0 ) ^ } ([C^ sin (0 + ip -<p ) |sin (0 + ~<t> )| 

- C 2 cos (0 + if#- 4> ) | cos ( 0 + ip _ <j>) | ]] ... (3.3) 

whore , C 2 referred to as the solar parameters in the analysis, 
are given by 

C 3 » l 3 A.(1 + P.) B .3 S/CmXC* ) 5 j = 1,2 ... (3.4) 

3.3. Effect of Solar Radiation Pressure on Attitude Dynamics . 

The governing nonlinear, non -autonomous differential 
Equation ( 3 . 3 ) does not have a closed-form solution, numerical 
techniques wore therefore used to evaluate the librational 
performance of the satellites in various practical situations. 
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Numerical integration was performed using Adams Extrapolation Method^ 
in conjunction with a suitable step-size of 2°. 

Figure 3 .2 shows the librational response in two typical 
situations 'with different orbital eccentricities. To bring out the 
influence of the SEP clearly, the response for = C 2 = 0, i.e., 
when the solar effect is ignored, has also been included for 
comparison. It may bo observed that the SRP increases the amplitude 
of librations when is positive. On the other hand when negative 
values are chosen for this parameter, the SRP seems to influence the 
system favourably reducing the maximum amplitude of librational 
motion. 

To understand the dynamical behaviour of satellites over a 
wide range of system parameters, numerous response plots were 
obtained. The resulting information is condensed in the form of 
system plots (Fig. 3*5 - 3 * 5 )* 

Figure 3.3 shows the effect of eccentricity on the librati- 
onal amplitude in throe different situations. It may be observed 
that in general the amplitude continually increases with increase in 
eccentricity. Hero, solar parameters , C 2 and the solar position 
angle 4 appear to be the other significant parameters of the system 
affecting the librational motion. 

Figure 3.4 indicates the effect of the solar parameters 
C 1 , C 2 on the librational response. It may be pointed out that in 
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caso of circular orbits and for = Cg = 0, the satellites with 
their long axes initially oriented along the local vertical do not 
execute the librational motion unless disturbed. In contrast the 
response results for any other non— zero combination of the solar 
parameters exhibit considerable amplitude of librations thus estab- 
lishing the detrimental influence of the SRP in the circular orbits. 
However, it is interesting to note that for elliptic orbits, the 
minimum amplitude does not correspond to the case = C 2 » 0, i.e., 
when the solar effect is ignored. Furthermore, through a judicious 
choice of the solar parameters, it is possible to reduce the libra- 
tional amplitude substantially. The significance of this result 
cannot be over-emphasized as it suggests an attractive possibility 
of countering the adverse influence of eccentricity on attitude motion 
of satellites. 

The effect of the varying solar position angle is indicated 
in Figure 3*5* It may be noted that the solar position has conside- 
rable influence on the librational amplitude of satellites, 
particularly in elliptic orbits. This suggests that a change in <j> 
signific.-mtly affects the choice of the solar parameters for achieving 
an improvement in the satellite librational performance. 

Thus, the numerical results presented here clearly establish 
the possibility of achieving considerable reduction of the maximum 
librational amplitude through an appropriate choice of the solar 
parameters. This beneficial aspect of the SRP appears to have been 
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either missed or entirely ignored so far. On the other hand, an 
impression seems to have gained ground that the uncontrolled solar 
effect is always detrimental to the attitude performance. This may 
perhaps partly explain why, in the recent developments of the various 
solar controller models, the emphasis has "been essentially on 
regulating the solar torque through the moving surfaces using complex 
feedback control systems. 

The system plots obtained here can no doubt be used for 
selecting the appropriate values of the solar parameters in several 
cases. However, to numerically generate the extensive data to enable 
the choice of those parameters covering all possible situations is 
neither feasible nor desirable. So, it was now proposed to treat the 
problem analytically with a view to develop the suitable criteria 
facilitating a judicious choice of the solar parameters. 

3.4 Synthesis of Control Strategy using Approximate Analytical Approach . 

It was felt that in view of the complexity introduced by the 

nonlinear and non-autonomous nature of the problem, it would be 

essential to make some reasonable assumptions for treating the 

equations of motion analytically. It is easy to see that for the 

practical situation of small amplitude librations assumed here, the 

16 

governing equation can be written as 

(1+ e cos ©W - 2eV sin 6 + 3^ sin^cos ip = 2e sin 6 + 

{ (i+e)^/(l+e cose)^} [G 1 sin( 8-<t>) |sin(© - ) I 

- C 2 cos ( 0 - <j> ) |cos( 6 - <P ) | 3 


... ( 3 . 5 ) 
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It may "bo pointed out that in general the deterioration of 

the attitude control performance characteristics is mainly due to the 

periodic excitation represented hy 2e sin 0 , Attempt was therefore 

made to explore the possibility of getting rid of this term with the 

help of the last two terms representing the solar pressure effects. 

To achieve this objective, the restructuring of the solar terms was 

considered so that with a suitable choice of the right hand 

side in the equation could be forced to disappear. Fortunately, 

Fourier expansion of the terms appearing within the square brackets 

in Equation (3*5) suggests that they can be effectively described by 

16 

their first harmonic . On incorporating this approximation, the 
Equation (5-. 3) simplifies to 

(1+ e cos 8)4> - 2ei/ sin0 + 3K^ sin^cos^= 2e sin 6+ (8/3 V 

{(l+e)/(l+o cose) } 5 (jC-| sin(9-4>) |sin(0-$)| 

_ C 2 cos ( 0 -<f> ) | cos (6 - <t> ) | II ... (3-6) 

On separately grouping the terms involving sin 6 and cos0 , it was 
found that the right hand side in the Equation (3*5) vanishes when 

C 1 = -0.75^ {e/(l+ e) 5 '} (1+e cose) 5 cos <J> 

C 2 = +0.75ir {e/(l+e) 5 } (l+o cos e ) 5 sin <}> ... (3*7) 

These results are of considerable significance, as they 
clearly bring out the effects of the orbital eccentricity as well as 
the sun— po sit ion angle on the desired instantaneous values of the 
solar parameters. It may be pointed out that by suitably regulating 
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the movement of the pairs of the solar surfaces and S 2 - S 2 

respectively along the z and x-axes (Pig. 3.1), & 2 aa d hence 

C 1f C 2 can he varied as desired. Thus using these relations, it 
may he possible to control the solar parameters suitably in order 
to overcome the adverse effect of eccentricity. 

A careful examination of the control relations reveals that 
the changing satellite position in orbit is likely to be of little 
consequence to the selection of , C 2 , particularly for low values 
of orbital eccentricities. Hence, these relations can be considerably 
simplified; 

C 1 = -0.75TT { e /(l+e) 5 } cos $ 

C 2 « +0.75TT (e/(l+e ) 5 } sin <f> ...(3.8) 

To examine the validity of the approximate analytical 
approach, the above expressions developed for , C 2 were substituted 
in (3«5) and the resulting equation of librational motion was 
integrated numerically. Some typical response results thus obtained 
are shown in Figure 3.6. It may be observed that the control schemes 
developed analytically are quite effective in limiting the amplitude 
of librations to much lower values. Even with the highly simplified 
control of solar parameters represented by (3-8) where 
C 2 = C 2 ( 4 >) , the amplitude reduction is found to be considerable. 

The more stringent control scheme represented by ( 3 - 7 ) * where 
C-j =C.j( 0 ,<|>), Cg = C 2 ( 0 » <{>) further improves the controller 


effectiveness . 
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To assess the influence of the important system parameters 
on controller effectiveness, the response results were analysed in 
details. Figure 3*7 shows the effect of the solar position angle on 
the system "behaviour. As the position of the sun changes continuously 
in the ecliptic plane, this study would he particularly useful in 
predicting the long range satellite lihrational dynamics. The analysis 
shows that the changing position of the sun does not significantly 
alter the attitude control characteristics, particularly for low 
values of eccentricities. 

Figure 3 *8 exhibits the effect of orbital eccentricity on 
maximum lihrational amplitude. It may be observed that both the 
control strategies are quite effective in reducing the amplitude of 
motion. For continuously increasing values of eccentricities, the 
effectiveness of the control schemes progressively declines. The 
degradation in performance is particularly severe with the simpler 
control scheme. This limits the applicability of the <j> -sensitive 
control to only slightly eccentric orbits. However, usefulness of 
the proposed, control scheme sensitive to both 0 as well as <j> , 
extends to much larger range of eccentricities. It is interesting to 
note that even for eccentricity as large as 0.3, it can limit the 
amplitude to well within 5°. l 11 view of the recently proposed use of 

highly elliptic geosynchronous orbits in order to maximize the net 

5 

on-station weight capability , this study may be particularly significant 
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Attention was now focused on the comparative assessment of 

the overall attitude performance characteristics associated with the 

two proposed control laws for the particular case of satellites in 

near-circular orbits. It was felt that, in this situation, even 

though the control operation sensitive to both 0 and 4 provides 

slightly better pointing accuracy, the <p -based action would result 

in considerably lower demand on the control system. It may be pointed 

out that the simple, 4 -sensitive control involving the single 

variable is easier to realize than the other one dependent on both 

0 and 4 . Furthermore, as the period of terns such as cos 4 and 

sin is one year -much higher than the orbital period of the 0-terms 

involved in the control relations, it is apparent that another major 

advantage of the simple 4 -sensitive control function lies in its 

much lower frequency. It may be possible to use this fact to 

advantage by effectively replacing the continuous periodic control 

action by a limited number of presjecified discrete operations over 

each cycle. To assess the effectiveness of this approach, the 

response plots were obtained for the varying degrees of discretization 

of the proposed 4-sensitive control criteria. The results of this 

analysis are shown in Figure 3.9. Here, 4 C represents the solar 

position angle with respect to which the solar control surfaces are 

\ 

assumed to be set initially and A <p denotes the maximum permissible 
deviation in 4 before permitting the resetting of the solar surfaces. 
It may be observed that the discretization leads to considerable 
deterioration in the attitude control characteristics. However, for 
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values of &<}> representing small changes in of say upto + 5° » the 
minor degradation in librational performance is likely to be •well 
within the tolerance limits for many space applications. This result 
is of considerable practical significance, as it establishes the 
feasibility of employing the solar surfaces suitably mounted on 
satellites to achieve substantial reduction in librational amplitude, 
even without exercising control during the long intervals extending 
upto as much as a ?/eek. The proposed control approach thus suggests a 
simple method for testing the concept of achieving librational 
control using the SRP with minimal changes in satellite design. 

3.5 Effect of Generalizing the Control Law on Satellite Attitude 
Performance . 

The open-loop control policies suggested here although quite 
effective in limiting the librational amplitude to low values, does 
not provide for damping the excited librations. It is, therefore, 
proposed to combine them with the proportional and derivative control. 
The resulting control law can be written as follows s 

c 3'°j1 + 0 32 " <b( 9 . •" b’® 

where 

= c 1 (4) or c 1 Ce , <f>) 

C 21 = C 2^ or ^2^' ’ ^ 

2k-nr<0 - $ <_ (2k+l) ^ 

+ (h 1 V+v if»); (2k +l) it < 9 + <t> <_(2k+2) ^ 
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+ (u 2 V+ v 2 ^) 5 (2k-t)iT< 0+\})_(J <_(2k+£) TT 

- ( U 2 ^ ' -4- v 2 4>) ; (2k+J)tr < <_(2k+-|> 


... ( 3 . 10 ) 

with | C j 2 | l.Cj max , k being an integer. 

The various results are obtained, by solving Equation (3 *9) > 
using the earlier numerical techniques in conduction with a suitable 
step size of O. 5 0 . 

Figure 3*10 shows the librational response of satellite in 
some typical situations. It may be noted that using the control relation 
Cj = Ck( ) obtained by the combination of ( 3 . 8 ) and ( 3 . 10), leads 

to significant reduction in the steady state librational amplitude; 
while (3.9) reduces the librational amplitude even further. 


To understand the behaviour of system and control parameters 
on controller effectiveness, numerous response plots were obtained. 
The resulting information is condensed in the form of system plots 
(Fig. 3.11 - 3*14). 


Figure 3 . 11 shows the effect of mass distribution on the steady 

state librational amplitude in several situations. It may be observed 

that the inertia parameter, has considerable influence on the 

amplitude for P. = v = 0 , i.e., when only open-loop control 
%} J 

C - = C • (<i> ) or C. = C.(6,4>) is employed. For the feedback control 

0 t) J J 

represented by nonzero values of y . and v., however, the change in K. 

0 J 3- 
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appears to have no effect on the system plots. Furthermore introduc- 
ing the additional terms to offset the effect of excitation induced 
hy eccentricity reduces the lihrational amplitude considerably. 

In general, the amplitude of oscillations is found to be lowest when 

the more stringent control strategy represented by C. = C . (0,<P , 4 ) 

«] j 

is employed. 

The influence of orbital eccentricity on the amplitude of 
motion is shown in Fig. 3.12. It may be noted that in general, 
increase in the eccentricity adversely affects the librational 
performance . However, the suggested modifications in the control 
strategy substantially improve the controller effectiveness and limit 
the amplitude to much lower values even in case of relatively large 
eccentrities. 


Figure 3* 13 shows the effect of control gains p. and v. on 
librational amplitude. It may be observed that a judicious choice 
of the gains can minimize the attitude variations considerably. Even 
for the best choice of controller gains, the suggested control 
modifications appear to bring about a considerable reduction in the 
amplitude of librational motion. It is apparent that even for 
eccentricities as large as 0 . 1 , the proposed controller can limit the 


amplitude approximately to within 0 . 1 °. 


Figure 3.14 shows the effect of 0. and the initial 

impulsive disturbance ip' * I n general an increase in C . affects 

o j max 

the attitude performance favourably to a point beyond which the 
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dependence "becomes rather weak. It may he observed that with the 
introduction of the suggested modification in the control law, it is 

possible to use much smaller values of C . 

Cmax 

The proposed controller presents an interesting possibility 
of controlling satellite orientation in elliptic orbits. The reali- 
zation of the control is easy and the size of the plates required quite 
modest. A preliminary estimate using the data of INTELSAT IV series 
of satellites shows that the plate size of 0.4 m with permissible move- 
ment of 20 cm. is adequate for orbital eccentricities of 0.1. For 
larger eccentricities, this requirement increases in nearly the same 
proportion. The semi-passive character of the system promises increased 
life-time with overall reduction in the cost. 

3*6. Concluding Remarks . 

The important features of the analysis and the conclusion 
based on them may be summarized as follows 

( i ) Contrary to as is generally believed, the SRP need not be 

always detrimental to the satellite attitude response. In fact, 
there are several situations, where for a proper choice of the 
solar parameters, the radiation pressure can substantially 
reduce the maximum librational amplitude of satellites in 
elliptic orbits. 

(ii) The system plots, developed numerically, summarize the satellite 
librational response characteristics over a wide range of 
parameters. Through these plots, it is possible to select the 
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appropriate values of the solar parameters minimizing the 
lihrational amplitude in several oases. 

(iii) The approximate analytical approach adopted here seeks to make 
use of the solar terms in the equation of lihrational motion 
to eliminate the major excitation effect of eccentricity 
normally responsible for the worsening of the attitude control 
characteristics of the conventional methods. The two control 
strategies thus evolved are found to be quite effective in 
restraining the satellite oscillations. 

(iv) The proposed control schemes are much simpler as they are 
dependent only on 9 and ip , which in turn can be expressed as 
a function of time alone. This makes them particularly suited 
for preprogramming of the control function unlike the earlier 
complex feedback control approaches sensitive to the lihrational 
disturbances and requiring on-board sensing and computation 

of the error signal for actuating the controls. 

(v) Using the relatively more stringent control law governed by 
the instantaneous values of the position of the sun (<J>) and 
satellite (®) leads to drastic reduction in maximum lihrational 
amplitude. It considerably extends the effective range of 
possible eccentricities for stable satellite operation. It is 
interesting to find that even for eccentricity as large as 0.5, 
the proposed controller can limit the maximum amplitude to 
well within 5°. This suggests the possibility of using it as 
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an auxiliary control device for improving the overall lihrational 
performance characteristics of satellites in highly elliptic orbits* 

The proposed control mechanism may thus have far reaching implications. 
In vi cy/ of its ability to overcome the adverse effect of eccentricity 
to a large extent ; it may lead to greater acceptance of highly 
elliptic orbits, particularly for multi-space missions. 

(vi) The simpler <J> -sensitive control strategy is quite effective for 

satellites in slightly eccentric orbits. Furthermore, it is possible 
to effectively replace this continuous periodic control action by a 
limited number of prespecified discrete operations ovor en.ch cycle. 

The analysis establishes the feasibility of employing the suitably 
mounted solar surfaces on satellites to improve the librational 
performance even without exercising control for long intervals. The 
proposed control approach thus suggests a simple method for testing 
the concept of librational control using the SRP with minimal changes 
in satellite design. 

(vii) Combining the suggested open-loop control with proportional and 

derivative feedback control appears to bring abmt a considerable 
reduction in the amplitude of librational motion. For example* the use 
of the modified control law enables the solar controller to limit the 

O 

librational amplitude to within approximately 0.1 , even for eccentri- 
cities as large as 0.1. 



4* ATTITUIE CQmOL OF SATELLITES IN NON-ST11IONAR.Y 2A-HOUR OB3ITS 


4*1 Introduction 

The advantages of having communications satellites in a 
perfect geosynchronous orbit are well understood. These orbits not 
only enable the use of simple ground stations but also facilitate 
the use of relatively simpler subsystems for station-keeping, 
attitude control and communications. It is, therefore, natural to 
conclude that any major deviations from the well established practice 
of using perfect synchronization between the satellite and the ground 
station is going to pose a major challenge on several fronts and a 
major effort is called for to meet these challenges successfully. 

In particular, the problem of using elliptic, inclined, 
24-hour orbits for communications is significantly complicated by 
continual periodic drift of satellites from their synchronous position 
directly over the ground station. The non-stationary nature of the 
satellite causing periodic oscillation of the line of sight about the 
local vertical leads to a continual periodic change in its orientation 
as viewed from the ground station. 

This problem can perhaps be easily handled through the 
active attitude control method. However, on-board energy requirements 
involved here are likely to offset the advantages of economy in 
booster capacity needed for placing the payload in the desired orbit. 
It seems, however, possible to overcome this difficulty through the 
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passive control schemes using solar radiation pressure. It is, 
therefore, proposed to explore the possibility of evolving a near- 
passive solar attitude controller which can vary the satellite 
orientation continuously so as to compensate for the steady or 
periodic satellite drift relative to the ground station. 

4*2 Formulation pf t.be Problem : 

Geometry of the satellites assumed to be moving in 
elliptic equatorial orbit is shown in Fig. 3 . 1 .. The angle, £ made 
by the z-axis with the line of sight defines satellite orientation 
as apparent from the ground station. Recognizing the relation 
between this apparent angle, £ and librational angle, ip as 

4 = £ + a ... (4.1) 

where 

a = tan -1 { sin(e - oj e t)/(r n - cos(e- <0^)) } ... (4*2) 

Using (4.1), the Equation ( 3 * 3 ) can be written as 
(l+e cos )£*' - 2e(l+£* ) sin + 3Fh sin Ceos £ 

= { (l+e)^/(l+e cos 9 )^ }[C^ sin (©+£ + a — <J> ) | sin(e+ £ + a - 4 ) | 

- C 2 cos(e+ £ + a - 4 ) I cos(6 + £ + a -4 )| 0 + e cose)a" 

+ 2e sin© a' + 3K. sin £cos £ - 3K i cos ip sin 4 ... ( 4 * 3 ) 

4.3 Effect of SRP on Attitude Dynamics ; 

To understand the dynamical behaviour of satellites over a 
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wide range of system parameters, various response plots were obtained 
numerically. The resulting information is condensed in the form of 
system plots (Fig* 4.1 - 4.3). 

Figure 4*"1 shows the effect of eccentricity on the 
librational amplitude in three different situations. It may be noted 
that the maximum librational amplitude increases with increasing 
value of eccentricity. The effect of solar parameters and C 2 on 
the librational motion is shown in Fig. 4*2. As before, it may be 
pointed out that the minimum amplitude does not correspond to the 
case when = 0. Furthermore, through a judicious choice of 

solar parameters, the librational amplitude can be reduced considerablj 

Figure 4*3 shows the effect of variation of solar 
position angle, 4> . It may be noted that <}> has considerable 
influence on librational motion. This shows that a change in cf> 
significantly affects the choice of the solar parameters to achieve 
improvement in the satellite librational performance . 

In view of the obvious limitations of numerical approach 
leading to a judicious choice of the various parameters, attempt 
was now made to use an approximate analytical approach for predicting 
the suitable choice of control parameters. 

4 .4 Analytical Approach % \ 


Due to complexity and nonlinearity of the nature of the 
problem several simplifying assumptions have "been made to treat the 
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problem analytically. For small arnnlitu.de motion, it is possible to 
approximate (-4.2) by the relation 

a Jfc sin (6- <o e t) (l + e oosfl)^ ... (4.4) 

where 

% = a n (l - e 2 ) - 1 

Linearizing (4 . 3 ) in conduction with (4 .4) , equation of 
motion can be written as 


n f 

(l+e cos 0 )£ - 2 e sin 0 £ + 3 K^ sin 5 cos £ = 2 e(l- 3 K^/a^) sin 6 

+ { (l+e)/ (l+e cos 6 ) } ^ [jj sin (9 - ( 4 > — ot )) 

| sin( 0 -( <p - a )) |- C2 cos(0 -(4 -a )) | cos(0 -(<}>- a )) | /] 

... (4.5) 

Following the steps as outlined in Chapter 5 > it can be 
shown that the effect of excitation due to eccentricity and gravity 
gradient moment can be counteracted by controlling the solar parameters 
according to the relations 


Cl = -0.75" c (1-3%/a.,) cos Y (l+e 00s e) 5 /0 + e) 3 

T Z 

C 2 = +O .75 ire (1-3%/a.j) siny (l+e c os9 ) ° 


... (4*6) 


where y = 4 - a anda % (0 -a^t)/a^ 


These results are of considerable significance, as they 
clearly bring out the effects of the orbital eccentricity inertia 
parameter and the sun position angle on the desired instantaneous value 
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of the solar parameters. It may be noted that by suitably moving 
the solar plates, the solar parameters can be varied in order 

to overcome the adverse effect of eccentricity. 


A careful examination of the control relations reveals 
that the changing satellite position in orbit is likely to be of 
little, cons&quence to the selection of particularly for low 

values of orbital eccentricities. Hence, these relations can be 
considerably simplified; 


C 1 = -0.75* o(l- cos4>/(l+e) 3 

= +0.75 * ©(l-jjKj/a.f) sin <J> /(l+e)^ 


... (4.7) 


To examine the validity of the approximate analysis, 
the above expressions for C^, Cg were substituted in (4*5) an h the 
resulting equation was solved numerically. The results are shown 
in Fig. 4 . 4 . It may be noted that the control schemes developed are 
indeed quite effective. However, the approximate control of solar 
parameters according to Cj = represented by (4*7) has almost 

the same effect as that due to Ch = denoted by ( 4 . 6 ). 

To understand the dynamic behaviour of system parameters 
on the controller effectiveness, the response results were analyzed. 
Figure 4.5 shows the effect of the solar position angle on the maximum 
librational amplitude. This study may be useful in predicting the 
long range satellite dynamical behaviour. The analysis shows that 
there is not much change in the librational amplitude especially for 
low eccentricities 4 
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Figure 4*6 shews the effect of eccentricity on the 

maximum libmtional amplitude. It may he observed that both control 

strategies are equally effective in reducing the amplitude of motion 

for the low values of eccentricities. As the value of eccentricity 

increases, the effectiveness of control schemes decreases. However, 

it is worthwhile to note that the amplitude can be limited well 
0 

within 5 even when e = 0.10. 

4*5* 'Effect of Generalizing the Control Law on Satellite Attitude 
Performance 

As the open loop control policies suggested above are 
not fully effective for higher values of eccentricities, it is 
proposed to combine them with the proportional and derivative control 
law. The resulting control law can be written as % 

Oj = C 51 + C j2 = C.(e,$, V) ... (4.8) 

where 

C in = C^,K ± ) or 0^0,*, K i ) 

C 21 = or O 2 (0,+, K ± ) 

C 12 = and C 22 = C 2 0f>,V) are given by (j.10) 

Since the modified governing equation of motion does not 
have closed-form solution, numerical techniques were applied to study 
the system behaviour. Hero, a stepsize of 0.5° was found to be suitable 
for numerical integration. 
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To examine the suitability of control given by ( 4 . 8 ), 

some typical response results were obtained (Fig. 4 . 7 ) . It may be 

observed that the solar control given by C . = C . (4 , 4* ) is quite 

*3 1 

effective in limiting the librational amplitude .The more stringent 

control schemes represented by C. = C .(4,IC. , 4, 4 ' ) combination of 

J «3 3- 

( 4 * 7 ) and ( 3 .IO) and C. = C.(9, 4 ,K. ,4,4* ) combination of (4*6) 
and ( 3 .IO) further improves the controller effectiveness. 

To assess the effect of system parameters on the steady 
state librational amplitude, numerous response results were obtained 
which aro presented here in a rather condensed form (Fig. 4*8). 

It may be observed that with the proportional and derivative control 
the steady state librational amplitude is reduced considerably. 

Thu amplitude can be reduced further using the modified control laws 
C- = C.(<f>, K. ,4,4') and C. = C .(6, 4, K.,4, V) . Furthermore, the 

J <J **- d o ■ L 

results suggest that the simpler control strategy denoted by 

C. = C . (4,K. ,4, 4 ) is quite adequate and should hence be preferred. 

4 . 6 . Concluding Remarks % 

Tho important features of the analysis and the conclusion 
based on them may be summarized as follows; 

^ -p'j Tho proper choice of solar parameters can reduce the maximum 

librational amplitude substantially . 

(ii) Tho system plots, developed numerically, summarize the 

satellite librational response characteristics over a wide 
range of parameters. Through these plots, it is possible to 
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bo loot the appropriate values of the solar parameters 
minimizing the librational amplitude 1 in several cases* 

(iii) The two control strategies 0. = C.(<J>,K.) and C . = C.(6, <f> ,K.) 

o 0 a 5 j i 

arc found to "be quite suitable in overcoming the major 
excitation effect of eccentricity and inertia parameter. 

(iv) The simpler control strategy represented hy C* = C-C 1 ?,]!.) 

j 1 i 

can he employed quite effectively for limiting the 
lihrational amplitude without appreciably compromising 
with pointing accuracy. Tho resulting solar controller can 
also ho used as an auxiliary control device for improving 
the overall lihrational performance characteristics of 
satellites in elliptic orhits if higher degree of pointing 
accuracy is required. 

(v) Tho addition of derivative and proportional control makes 
tho controller more effective. This generalized control 
limits the steady state lihrational amplitude to well within 
5° for eccentricities even as largo as 0.20. 

(vi) Tho generalised control policy indicated hy ) 

appears to he tho host as it generally leads to minimum 
limit cycle amplitude of lihrations. 



5. CLOSING- EEMifflKS 


5.1 Summary of the Conclusions 

As already stated the study aimed at examining the various 
aspects of motion related to the use of 24-hour, non-stationary orbits 
for satellite communication. Tho important conclusions based on the 
analysis can bo summarized as follows”. 

(i) The use of the proposed non-synchronous orbits can enable 
substantial saving in the booster capability requirements. 

(ii) Even for relatively large eccentricities, and with moderate 
orbital inclinations and ground station latitudes, the 
uninterrupted visibility is possible. 

(iii) The requirements of peak power and total energy for 
communications show a continuous increase with increasing 
orbital eccentricity. However, this may not pose a major 
limitation, particularly with proper scheduling of 
communication periods and when this main objective 
satellite communication objective is combined with various 
other important missions, e.g., remote sensing for survey 
of earth resources and better weather -forecasting 
scientific space research etc. 

(iv) To provide for the changing satellite attitude so as to 

achieve its fixed orientation with respect to the moving 
line of sight appeared to be a major challenge involved 
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involved, in the use of the proposed nonsynchronous orbits 
for communications . The investigation showed that this 
challenge can he successfully net by employing the passive 
attitude control method using the solar radiation pressure. 

5*2. Recommendations for Future Work 

There are several aspects related to the use of 24-hour 
non- synchronous satellites for communications which have not been 
considered in the present investigation. It would be worthwhile to 
carefully examine those aspects which may be summarized as follows: 

(i) The ground station may require using more sophisticated 
system. It would, therefore, be essential to examine the 
effect of the resulting increase in the cost of ground 
component on overall economics. 

(ii) A transmission power control will become necessary to 
counteract the variation of the range caused by orbital 
ellipticity. 

(iii) Several communication problem related to frequency coordina- 
tion, intermodulation are likely to be adversely affected. 
Besides, it would be desirable to extend the analysis to study 
the attitude control characteristics of the proposed solar 
model as affected by the orbital inclination, roll and yaw 


motion 
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APPENDIX 


Calculation of total velocity increment in transfer of satellites 
from parking orbit to 24-hour coplanar elliptic orbit. 

The geometry of satellite motion and its tr ans fer from 
circular parking orbit to the final 24 -hour elliptic orbit is shown 
in Pig. 2.1. It is easy to show that the circular velocity at P 
is given by 

V - (a.i) 

c 

The corresponding velocity at P along Hohmann transfer 
trajectory is given by 



The velocity increment required at p for injecting the 
satellite from parking orbit into the 24 -hour elliptic orbit is 
therefore given by 


av g = 4 /ag ^ a 2 / a 1 H 

(A. 3) 

On non-dimensions Using the velocity increment, we get 

AV e = (A .4) 


AV. 


where 


A V 4/a 2 
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A .2 Calculation of total Telocity increment in transfer of satellites 

from parking orbit at an inclination with respect to the equatorial 
piano to geosynchronous orbit. 

As shown earliest, velocity increment required at P to change 
the satollitu orbit to Hohmann transfer trajectory can be written as 

Av p = ^ V /a 2 Q/2a 2 / ai - 2a 2 /(a 1+ a 2 ) - / a^/a"] (A. 5 ) 

Similarly, tho velocity increment required at A in order to inject 
the satellite into geosynchronous orbit can be obtained 

&V A " ^/ a 2 CAl- ( a i/ a 2 ) /2a 2 / a 1 " ^ a 2 / ( a -j + a 2 ) }2 +(2sin i/2> 2 

(a. 6) 

where i is inclination of parking orbit, combine (A. 5 ) and (A. 6) 
results in the following non-dimensional expression for the total 
velocity increment 



(A. 7) 

where 

AV C - AV C / Ai/ag 
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Fig. 2-1 Geometry of 24-hour circular and elliptic orbits 














Fig. 3-5 System plots showing the effect of 
solar position angle on maximum 
librational amplitude 
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Fig. 3 6 Typical satellite response showing the 

attitude control through the proposed 
solar controller for (b) e =0.1 



i max 


56 



0 120 240 360 


4 >° 

Fig, 3*7 System plots showing the attitude control 
characterstics as affected by the solar 
position angle. 
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Fig. 3*8 System plots showing attitude control 
characterstics as affected by the 
eccentricity- 
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Fig.3-9 System plots showing the effect of ovinyir 
solar position angle from its nominal value 
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Fig. 4-1 Librational behaviour of satellite 
relative to the ground station as 
affected by eccentricity 




ntional behaviour of satellite relative to tl 
ons as affected bv solar narnm(?te>rQ 




3 4-3 Librational behaviour of satellite 
relative to the ground as affected 
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Kj - 1*0 e ~ 0-1 4* ~ 1 0° 

Y o = 0 Y--0 


Cj = C; (4>,Ki) 
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Fig. 4-4 Typical satellite response showing the 
control of librations of satellite about 

line of sight as achieved through the 
proposed solar controller 



CjrCj (4>, Kj) 

Cj = Cj (0 ,<f> »Kj ) 


Fig. 4-5 Librational control of satellites relative 
to the ground station as affected by 
position of the sun 



Fig.4-6 Librational control of satellite relative 
to the ground station as affected by 
the eccentricity 
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Fig. 47 Typical response plots showing the relative 
attitude control performance 
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Fig. 4-8 System plots showing the effect of eccentricity and 
control gain jj on attitude control performance relative 
to the ground station 
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*IBFTC RKADE 

SUBROUTINE RK ALL C X, XL , Y , Y 1 , U 


(. 

r 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


*»«•«#»*****» tt* ##«**«*«««#»**«*«#' ; a 


* SUBROUTINE F OK SOLVING THE SECOND ORDER « 

* ORDINARY DIFFERENTIAL EGl AT ION-- * 

* Y» * = F( X, Y, Y* ) 

* X--AN ARRAY Oh X-AXIS VALUES « 

* XL — FINAL VALUE OF X * 

* Y — AN ARRAY OF Y-AXIS VALUES 

* Y1--DY/DX * 

* H1--STEF LENGTH * 

* II— NO. OF LINES ALREADY PRINTED AND « 

* LEFT AS BLANK « 

* JW— INPUT-OUTPUT DEVICE FCR PRINTING « 

* KK=10 IF PRlNflNG OF VALLES IS NOT « 

* REQUIRED « 

* NN--FLAG INDICATOR * 

* =0 IF X AND Y IN RADIANS * 

* =1 IF X AND Y IN DEGREES * 

« =2 IF X AND Y IN RADIANS AND INTEGR- * 

* AT ION IS F RON G TC 2*PI(0R MULTIPLE * 

* OF PI ) * 

* FUNCTION F TO BE PROVIDED BY USER « 

* ADAMS EXTRAPOLATION METPCC USED TO SCLVE « 

* DIFFERENTIAL EQUATION « 

* RUNGE KUTTA FOURTH ORDER METHOD USED FOR # 

* STARTING VALUES * 




REAL K1»K2»K3»KA 

DIMENSION FX(6> ,DfcLl<6),EEL2<6>,DEL3*6) ,EEl.4(6> ,CEL5(6) 
DIMENSION X(L)»Y(L),Y1(L) 

COMMON/ VKJ/H1, 1 t*oW,KK,NN 
P I = 4 .*ATAN(1 . ) 

PP=PI/180. 

PI2=2.*PI 
H = Hi 

IF (NN.EQ.l) H=PP*Hl 
H2~H*H/2 . 

H3=H/2. 

N =5 
M = N+1 
C 

C **« RUNGE-KUTTA ME I HOD *** 

C * 

IF (NN.EQ.l) Y=PP»Y 
DO 1 1 = 1, N 

IF(NN-EQ.l) X(i)-PP*X(I) 

FX(I)=F(X(I>,YU),YKI>) 

K1 = H2»F(X<I),Y(1),YKI)> 

K2 = H2#F < X C I ) +H3 > Y U ) + f-3*Yl(I)+Kl/4. » Y l { I ) + K1/H > 

K3 = H2#F(X(I)+H'5.Y( I )+H3«Yl U5+K1/A, ,Y1.( I5 + K2/H) 

K4=H2*F (X( I ) + H, Y( 1 ) +P«Y1( I )+K3» Y1 ( I I+K3/H3 ) 

IF(NN.EQ.l) X ( 1 ) = X ( l ) /F F 
K - I + 1 

Y(K)*Y( I)+H«Y1( i> +(Kl +K2+K3 >/3. 
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V 


YKK)=Y1( I) + t K1 +2 .«K2+2 .«K3 + K4) / 3 ./H 
3 X ( K ) = X ( 1 ) +H 1 

IF(NN.EQ.l) X(K) = PP*X(K) 

FX(K)sF(X(K), YCK) ,YKK> ) 

IFC.NN.EQ.1) X(K)=X(K)/FF 

C * 

c »** AC AYS EXTRAPni. A! ION YETKE « *« 

C 

2 DO 3 1=2, M 

3 DELHI )sFX(I )-FX(I-l ) 

DC 4 1=3, M 

4 DEL2C I > -DELI ( I ) -DELI ( 1-1 ) 

DC 5 1=4, Yl 

5 DEL3( I ) =DEL2( I)-DEL2< 1-3 ) 

DO 6 1=5, M 

6 DEL4( I ) =DEL3( I)-DfcL3( 1-1) 

DEL5(M)=DEl.4(M)”DfcL4 (K-l ) 

K = K + 1 

Y(K) = Y<K~l)+H*Yl(K-3.)+MP*(FX(H>/2.+nE!_1(Y)/6.+EFL5fM>/8,+ 
1DEL3(M)#19,/180 . + DEL4 (Y )«3,/32, +DEL5' •-') «P63 ./IOORO. ) 

Yl(K) = Yl(K-i)+H*(hX(Y )+CELl(M)/2.+DbL2(H)«5;/12.-‘DEL3(P)«3 ,/8 .+ 
lEEL4(M)*251./72U.+CEL5(Y)«95./288. ) 

CO 7 1=1, N 

7 FX( I ) =FX ( 1+1) 

X(K)=X(K-1) +H1 
IF(NN.EQ.l) GO 10 8 
IF(NN.EQ.2) GO 10 9 
FXCY)=F(X(K),YUO,YKK) ) 

IKX(K)-XL) 2,2,10 

8 XX=X(K) 

X(K)=X(K)-FLOAHIFIX(X(K)/360,))»3 60. 
ir(NN.EQ.l) X(K)=PP*X(H) 

FX(H)=F(X(K),Y(K),Y1(K)) 

X ( K ) =XX 

IF (X(K)-XL) 2,2,10 

9 XX=X(K) 

X(K)sX(K)-FL0AI(IUX(X<K>/PI2))*PI2 

FX(Y)=F(X(K),Y(K),Y1(K)) 

IF ( X ( K ) -XI. ) 2,2,10 

10 I F ( NN , NE . 3L ) GO 10 12 
DO 11 1 = 1, L 

11 Y( I )=Y(I )/PP 

12 IFCKK.EQ.10) RhIUHN 
C »•»« 

C *** PRINTING THE VALUES CF X AND Y *** 

C 

Jl=l 
J2 = 6 

13 11=11+1 

I F ( 1 1 , EG . 70 ) GO 10 16 
' WRHE(JW,155 (X(I),Y(I),I = ul,J2) 

*J1 = Jl+6 
J2= J2+6 

IF(J2-L) 12,14,14 

14 vl2 = J2”6 + l 
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WRITE(JW.i5) (XII), Y(I)# I*v2»L) 

15 FCRNAT(1X,6<HU.1,HQ.5>> 

RETURN 

16 WR i TE ( J W > 17 ) 

17 FORMAT ( 1H1 » // ) 

11=0 

GO TO 13 
ENT 


alBFTC MAX 

SUBROUTINE M AX C Y , N , YM A X ) 

C «•«■»««•»**«•»•»#*■»*■»•»« •s*******'® 1 *'*********' 8 ' 0 * ■s « 

C * SUBROUTINE USfcU TO FINE Tl-E MAXIMUM * 

C * ABSOLUTE LIBRA I I ON AL AMFLITUDE * 

C * Y--INPUT VECTOR * 

C * N - ~NO , OF ELEMENTS IN VECTOR Y * 

C * YMAX — MAXIMUM ABSOLUTE CF VECTOR Y * 

(3 ###&##*#####•&#****#*»**#* ***#### 

DIMENSION Y ( N ) 

DO 1 1=1, N 
Z=A8S ( Y ( I ) ) 

1 IF (Z.GT.YMAX) YMAX=Z 
RETURN 
END 


« IBFTC FS 

FUNCTION F ( THE ! A , S 1 , S 1 1 ) 

C ** ** 

C # FUNCTION SUB-RKOGRAM ISED TO SOLVE THE « 

C * GIVEN DIFFERENCIAL EOLATIONS USING * 

C * CONTROL POLICIES FOR SATELLITE TO EXECUTE » 

C * ABOUT THE LOCAL VERTICAL « 

C * THETA — POS I T 1 ON OF SATELLITE * 

C * si — ubrational ANGLE s 

C * S II — DERI VATIVE OF L IBR AT ICNAL ANGLE WITH * 

C » RESPECT TO ANGIE THETA * 

C * E — ECCENTR I C I i Y * 

C * KI — INERTIA PARAMETER * 

C * PHI — SUN POSIliON ANC-LE « 

C * P I --CONST ANT 13.14159265) « 

C * AMU1,ANU1,AMU2,ANU2-CCNTRCLLER GAINS « 

C # C1MAX, C2MAX--MAX IMUM PERMISSIBLE VALUE OF * 

C * SOLAR PARAMETERS * 

C * LL— FLAG INDICATOR TC ISE ONE CONTROL * 

C « STRTAGEY A 1 A ONE TIME * 

C * Cl, C2 — SOLAR PARAMETERS * 

REAL KI 

COMMON E, KI » PHI >Pl,AMUl,AMJl»AMU2,ANU£,ClMAX»C2MAX»LL 

PI2=2,*PI 

A=1 , +E 

B = 1 , + E*COS ( THE ! A ) 

C* E*S I N C THET A ) 
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C=THETA+S1-PHI+HI2 
D=D-FLOAT( IFIX(L/P12))«FI2 
IF(LL) 10,20,30 

10 ci=o. 

02 = 0 , 

GC TO 40 

20 Cl=-0,75*PI*t*cnS(PHI>/A««3 
C2= 0 ,75*PI#fc»S IN (PHI )/A««3 
GC TO 40 

30 Cl = -Q .75*PI*E#U0S(PHI )« <E/A >»*3 
C2 = 0 ,75*PI#E*S1N(PHI )« (E/A )**3 
40 C11=-CAMU1*SU + ANU1«S I) 

IF(Cll.GT.ClMAX) C11=C1N AX 
IF ( C11.LT. ( -C1M AX ) ) Clls-CIPAX 
IFCD.GT.PI . AND. L ,LE .2 ,*P I) C11=-C11 
C1=C1+C11 

C21 = (AHU2«Si 1+ANU2#SI ) 

I F ( C21 , GT , C2H AX ) C21 = C 2NAX 
IFIC21.LT, C-C2MAX)) C21=-C2NAX 
IF(D.GT.PI/2..ANU.C.LE.3.«PI/2.) C21=-C2t 
02=02+021 

G = (2,*C*(1.+SU >-3.*KI*SIMSI )*CQS(SI ')/B 

H=A«*3/B*»4*(C1«S1NCE>*( AES(SIN(D) ) ) -C 2*CC5 < 0 > * ( APS ( CCS ( C > ) ) ) 

S 12 =G + H 

F= S I 2 

RETURN 

END 


* IBF T C FZ 

FUNCTION F ( THE I A , S I » S 1 1 ) 

C « 4MH* ## ****** ** ** ** ** ** *# * # # # ** ** *45 * 

C * FUNCTION SUB-PROGRAM USED TO SOLVE THE * 

C * DIFFERENTIAL EQUATIONS USING DIFFERENT , * 

C * CONTROL POLICIES FOR SATELLITE TO EXECUTE * 

C * ABOUT THE LINE OF SIGHT c 

C * THETA— POSITION OF SATELLITE * 

0 * SI— librational angle * 

C * S II — DERIVATIVE OF LIERATICNAL ANGLE WITH « 

C * RESPECT 10 ANGLE THETA « 

C * fc— ECCENTRIC! I Y * 

C « KI — INERTIA PAKAMEIER • 

C * PHI --SUN POSIliON ANGLE * 

C * P I --CONST ANT 13.14155265) * 

C * AMU1»ANU1, A MU 2, ANU2— CONTROLLER GAINS * 

C * C1MAX.C2MAX — MXIPUP PERHISSIBLE VALLE OF * 

C « SOLAR FARAHETERS * 

C * AN — CONSTANT ( 6 . 628 ) « 

C * LL--FLAG INDICATOR TO USE CNE CONTROL * 

C « STRTAGEY A! A ONE TINE * 

C * Cl, C2— SOLAR PARAMETERS « 

C « Wl — THETA-WE# l * 

C * WE--ANGUL AR VELOCITY CF EARTH * 

C * T — TIME * 

C **#**********4f^***^****4?4l«|l«4i*lf****^*^*^45*****^ 
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REAL KI 

COMMON E » K I » PH ; * H I , AH 1 , A NL 1 , AMU2 , AM 2 , C1M A X , 02 M A X , AN . L L 
PI?=2,*PI 
A = 1 . +E 

8 = 1 . +E*CQS ( THE I A ) 

C= E*S IN ( THE FA ) 

d=theta+si-phi 

D = C-FLOAT( IF IX (C/P 12) )«FI2 
G = 1 , »E*E 

ai=an#g-i, 

IF < LL > 10*20,30 
10 C1=Q. 

02 = 0 . 

GO TO 40 

20 Cl=(-0,75*PI*t*COS(PP I ) + 2 .25«PI*KI/A1»E*CCS(PHI ) 

C2= ( 0 .75*PI*E»SIN<PH I )-2 . 25»PI «K I /A1*E«S IN (PHI ) ' 'A«*3 
GO TO 40 
30 F=1 ,-E 

X = ( SORT ( A ) -SQR I <F) ) *S I N ( T P ET A ) / ( 2 . * ( SORT C A >*CQS ( THE T A/2 . )«*2+ 
1SQRT(F )*SIN (THE 1 A/?. )««2 > > 

Y=ARSIN(2.*X/(1.+X«*2) ) 

W1=SQRT(G)/8*C+Y 
H^PH I -Wl/Al 

Cl=(-0.75*PI*E«COS(H)*2.25«PI*KI/Al*E*C0S(h ) >«(B/AT«*3 
02= ( 0.75*PI*E*3IN(H)-2 . 2 5#P I »K I / A1 *E «S I N C H ) )*(B/A)**3 
40 CU=-<AHU1*SI1+ANU1«SI) 

IF(CU.GT.CIMaX) C11=C1MX 
IF(CU.LT.(-C1MAX)J C11=-C1PAX 
IKC,GT.PI.AND.U.Lfc.2.«Fl) C11 = -C11 
C1=C1+C11 

021= (AMU2*SI1 + ANU2*SI ) 

1F(C21.GT.C2MAX) C21=C2MAX 
IF (021 .LT. (-C2MAX) ) C2l=-C2MAX 
IF(C.GT.Pl/2..ANC.D.LE.3.«FI/2.) C21=-C2l 

02 = 02+021 

G = (2,*C*(1.+SU )-3.*KI«$IMSI >*C0S(S1 t ) /B 

H = A*»3/8**4*(C1*SIN (E >* < AES <SIN<D) ) ) -C2*C0S < D ) * ( ABS< COS (D > ) ) ) 

S I2=G + H 

F=S 12 

return 

END 


« IBFTC ZIVAL 

SUBROUTINE Z I V AL ( T HET A 1 , 2 I , N ) 


C » SUBROUTINE USED 10 EVALUATE THE VALUE OF « 
C # ANGLE ZI FROM LiBRATICNAL ANGLE SI « 
0 * THETA1 — PQSIT1UN OF SATELLITE * 
C * ZI— INPUT VECrOR SI AND RETURNS AS OUTPUT « 
C # VECTOR ZI * 
C * ALPHA--ANGLE BETWEEN TEE LOCAL VERTICAL * 
0 • AND THE LINE CF SIGHT * 
0 • N — NO. OF ELEMENTS IN VECTOR ZI * 


C #####*<**# #***######**###****##***#**#'##«i#**tt#**# 



REAL K I 

DIMENSION THbTU(N5,2I (N) 

COMMON E» K I .PNC »FI , AM 11, AMI, AMU2 , ANU2 , C1MAX , C2MAX- AN 

P 1 2 = 2 . *P I 

A=1 . +E 

F=1 .-E 

6=1 . -E*E 

DO 10 I =1 • N 

THETA = THETA1U ) 

B=1,+E*C0S(THEIA> 

C= fc*S I N ( THET A ) 

X = < SORT (A )-SQPI (F ) ) *S IN ( T HE T A ) / ( 2 . * (S GR T < A ) «CQS (THETA/2 .)**2+ 
!SGRT(F)«SlN(THt ! A/2. )**2 ) ) 

Y = A RS IN (2 ,*X/ (1 . + X* #2 )) + F 12 

W1=SQRT(G)/B*C + Y 

W1=W1-FL0 AT ( I FI X( Wl/F 12 ) ) *F 12 

rn=an*g/b 

ALPHA = ATAN(SIN( W1 )/ (RN-CCS(K1 ) ) ) 

10 2 I ( I ) =Z I ( I > -ALPHA 
RETURN 
END 



